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In the paper we present the qualitative analysis of rolling motion without slipping of a homogeneous round disk on a 
horisontal plane. The problem was studied by S. A. Chaplygin, P. Appel and D.Korteweg who showed its integrability. 
The behavior of the point of contact on a plane is investigated and conditions under which its trajectory is unit are 
obtained. The bifurcation diagrams are constructed. 




1. Introduction 

For the first time the motion of a heavy dynamically symmetrical round disk on a horizontal absolutely 
rough plane was investigated by GSlesser (1861) |28| . N.Ferrers (1872) [5], K.Neumann (1886), and 
A. Firkandt (1892). These studies eventually (after unsuccessful attempts by Neumann and Lindelof) 
lead to the correct form of equations of motion. This form differs from the usual (Lagrangian or 
Hamiltonian) equations of mechanics because of the nonholonomic constrain showing that the velocity 
of the point of contact of a disk with a plane is zero. We shall not discuss in detail the general forms of 
the equations of the nonholonomic mechanics (they are presented, for example, in [2Sj, [IS])- Instead, 
we concentrate on the pretty obvious form of these equations obtained from the general principle of 
dynamics — the conservation law of the moment of momentum written in the disk- fixed axes. 

S. A. Chaplygin (1897) was the first to show the integrability of the problem on rolling motion of 
a disk. He presented the reduction of the problem to the analysis of hypergeometric quadratures in 
paper [Hj, where he showed also the integrability of the problem on rolling motion of an arbitrary heavy 
dynamically symmetric body of rotation on a horizontal plane — in the latter case the problem is 
reduced to the integration of the linear differential second-order equation. The integration of equations 
of motion of a disk in hyperelliptic functions was also performed in 1900 independently from each other 
and from Chaplygin by P. Appel [2] and D. Korteweg |14j . Sometimes the problem on rolling motion 
of a disk is referred to as Appel-Korteweg problem (or simply Appel problem), but this is, probably, 
not quite correct. In 1903 the same result has been rediscovered by E. Gellop however he used 
the Legendre functions. 

Despite of the explicit hypergeometric quadratures the various qualitative properties of disk mo- 
tion were not studied for the long time. There were mainly studies of stationary motions and of 
their stability (the corresponding bibliography is presented in book |2B])- Some qualitative properties 
of the disk motion have been discussed only in papers S. N. Kolesnikov JHJ and Yu. N. Fedorov 
The first paper shows that for almost all initial conditions the disk never falls onto a plane and 
the second one present the procedure of investigation of the reduced system. Analogous results 
for the dynamically asymmetrical disk and disk moving on an inclined plane (nonintegrable prob- 
lems) were obtained in pQ, ^Sl- Among the modern works analyzing the rolling motion of the 
disk we shall note papers O.M. O'Reily [221, R- Cushman, J.Hermans, D.Kemppainen [Jj, and 
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A. S. Kuleshov [U devoted to the study of bifurcations and stability of stationary motions of the 
disk. 

General results of a qualitative analysis for the rolling motion of a heavy body of rotation were 
obtained in paper N. K.Moshuk [21] ■ The paper include the frequency analysis, application of the 
KAM-theory, and basic qualitative properties for the motion of the point of contact. It appears that 
the point of contact performs the composite bounded motion: it periodically traces some closed curve 
which rotates as a rigid body with some constant angular velocity about the fixed point. Thus the 
realization of some resonance relation between frequencies makes possible the drift of the body of 
rotation to the infinity. 

In this paper we develop these qualitative considerations and complement them with the 
computer analysis. We also present various types of trajectories which are traced by the point 
of contact in the body-fixed and relative frames of references since they have curious forms 
which are difficult to predict. Using the computer modelling we explicitly investigate the hy- 
pothesis about the drift to the infinity under the resonance conditions. We present the most 
general three-dimensional bifurcation diagram in the space of the first integrals and the com- 
plete atlas of its sections by various planes, constructed with the help of computer mod- 
elling. 

In this paper we also present a new method of reduction of the problem to an one-degree integrable 
Hamiltonian system and explicitly consider the existence of Hamiltonian formulation for different 
variants of equations of motion of the problem. 



2. The rolling motion of a rigid body on a plane 
2.1. Equations of motion and their integrals 

Let the rigid body in an exterior field of force perform a rolling motion on a plane without sliding. In 
this case the equations of motion have the most convenient form in the body-fixed frame of references 
which axes are directed along the principal axes of inertia of the body and the origin is situated at 
the center of mass. In the following text all vectors are assumed to be projected on these axes. 

The condition of absence of slipping thus becomes 

v + usxr = 0, (2.1) 

where v, u are the velocity of the center of mass and the angular 
velocity of the body and r is the vector directed from the center of 
mass to the point of contact (see fig. . 

Let's denote the projections of the fixed basis vectors to the 
moving axes by a, (3, 7 (the vector 7 is perpendicular to the plane) 
and by (x, y) we shall denote the coordinates of the projection of the 
Pig 1 center of mass onto the plane in the fixed frame of references. We 

assume that the field of force is potential with a potential depending 
only on the orientation of the body U = U(a, (3, 7). The complete set of the equations of motion 
defining the given system can be represented in the form 

M = Mxuj + mrx(uxr) + ax^- + (3x^- + jx (2.2) 

da. 0(3 07 

a = a x uj, (3 = (3 x u, 7 = 7x0;. (2.3) 

The equation ()2.2j) describes the evolution of the vector of moment of momentum for the body with 
respect to the point of contact M and (|2.MJI concerns the evolution of the fixed basis vectors in the 
body-fixed frame of references. 
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The motion of the center of mass can be obtained in quadratures from solutions of the equa- 
tions dH2), (|2T5)) as follows 

i = (rxw, a), y = (r x u, (3). (2.4) 

The expression of the vector of moment of momentum with respect to the point of contact M 
can be written in the following form 

M = Iuj + mr x (u x r), (2.5) 

where I = diag(ii, I2, I3) is the tensor of inertia of the body. In turn r can be uniquely expressed 
(for a convex body) through the normal to the plane 7 from the equation 

7 |VF(r)|' 

Here F(r) = is the equation of the body's surface. 

Let's consider a motion of the point of contact on a plane. If we denote the position of the point 
of contact on the plane in the fixed frame of references as (X, Y) , then the equation of motion for the 
point of contact can be presented in the form 

X = (r,cx), Y=(r,f3). (2.7) 

where r is determined from equations Q2.2JI - (|2.6|) . Actually X and Y are projections of the velocity 
of the point of contact in the relative frame of reference onto the fixed axes. 

The equations of motion in the form similar to Q2.2|) - ()2.3|) are presented, for example, in book jllj . 
They can be obtained also by means of Poincare - Chetaev formalism [3j with undetermined Lagrangian 
coefficients; these coefficients shall be eliminated with the help of the constrains' equations (|2.1|) . 

The system (|2.2|) - Q2.3J) generally has seven independent integrals of motion, six of them are 
trivial geometrical integrals: 

« 2 = 1, /3 2 = 1, 7 2 = 1, (2 8) 

(a, 0) = 0, 03, 7) = 0, (7, a) = 0. 

The seventh is the integral of energy 

i(M, w) + U{a, (3, j) = h = const. (2.9) 

Generally the given system has no other additional integrals and the possibility of its integrability in 
concrete cases depends on the presence of additional tensor invariants (measure, fields of symmetry, 
integrals). 

2.2. The rolling motion of a heavy disk 

Let's consider the case of rolling motion for an axially symmetric disk of radius R in the field of gravity. 
The field is, obviously, also axially symmetric with the potential depending only on 7. Moreover, we 
suppose that the disk is dynamically symmetric, i. e. l\ = Ii- The potential energy in this case has 
the following form 



U = -mg(r, 7) = mgR^l (2.10) 

The equation of surface for the disk is F(r) = r\ + r% — R 2 - Substituting it in the equation (|2.fi|) 
and solving with respect to r we obtain 

n = . r 2 = , r 3 = 0. 2.11 
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As the potential energy depends only on 7, in the equations of motion ()2.2j) - ()2.3|) we get the 
separate system of six equations. 

M = M x u + mr x (u x r) + mar x 7, , 

v ; ' (2.12) 

7 = 7x11), 

Expressing tj, r from relations (|2.5|) . (|2.11|) we shall get the closed system for the variables M, 7 
similar in many aspects is to the Euler-Poisson system in the Lagrange case, however the obtained 
system is much more complicated than the last one. 

The equations (|2.12|) preserve the geometrical integral 7 2 and the energy (|2.9j) . in addition they 
allows the standard invariant measure (with a constant density). For the integrability (by Euler- 
Jacobi ^Jj) of these equations we need two additional integrals. In the following we describe the 
method of derivation of these integrals. 



The possibility of separation of the system (|2.12j) from the general system (|2.2[) - (|2.3jl is connected 
to the symmetry with respect to the rotations about the vertical axis defined by the vector 7. The 
system ([2.12)1 is invariant with respect to the field of symmetries commuting with the vector field of 
the problem. 

* = ai A - + 012 A - 02 A + a3 A - ^ (2 - l3) 

It is possible to show that the variables M, 7 are the integrals of field (|2.13|) that is v^(Mi) = 0, 
%(7i) — 0, % = 1, 2, 3. According to the general Lie theory variables M, 7 define the reduced 
system. For the classical Euler-Poisson equations the corresponding reduction is the Raus reduction 
with respect to the cyclical angle of precession. 

In addition to the field of symmetries (|2.13j) the equations of motion (|2.2|) - (|2.3|) for the axially 
symmetric body allow one more field of symmetries corresponding to the rotation about the axis of 
symmetry of the disk. 

V in = M^ — Mo— h Ti — T9— — h 

(2.14) 

+ai £- 2 ~ a2 A +l3i w2~ P2 Wi- 

It is possible to show that integrals of the field (|2.14j) are projections of the moment and normal 
to the plane of disk onto the fixed axes of coordinates 

N = ((M, a), (M, 13), (M, 7)), n = (a 3 , /3 3j 7s)- 

The equations of motion for these variables can be presented in the following form 

N = mr x (u> x r) + mgf x n, ^ 
h = u} x n, 

where symbols u>, r denote the same vectors, but projected onto the fixed axes (that is Zj\ = 
= (oj, a), . . . , r\ = (r, a), . . .). The explicit expression of the components of the vector r is 

/ Ra 3 -/ 3 -R/?373 



Rjl-li . (2.16) 



The vector N is expressed through u by the formula 

N = I\u> + (I3 - h)(u>, n)n + mr x (r x r). (2.17) 
Remark 1 . Such reduction is also possible for an arbitrary body of rotation. 
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2.3. A reduction to the integrable one-degree Hamiltonian system 

Let's describe the process of reduction of order with respect to the both fields of symmetries (|2. 13|) and 
(|2.14|) . For that we shall choose the simultaneous integrals of these fields as variables of the reduced 
system. According to jH], the most convenient algebraic set of such variables is 

73, K\ = M l7 i + M 2l2 =N 3 - 73 (JV, n), 



K 3 = 71 M 2 - 72 Mi = Nin 2 - N 2 n\. 



The equations of motion in the new variables become 

K 3 



73 



h + mR 



2 ' 



K\ = 13 — K 3 K 2: 

(h + mR 2 )y/h{I 3 + mR 2 ) 

ir mR 2 K 3 K X 



(h + mR?)^/h{I 3 + mR?) 1 - 73 ' (2 ' 19) 
1 ~ 7,1 \ 7l J i + mR 



^h{I 3 + mR 2 ) 



K X K 2 + mgR^ J 1 - t| 



The equations (12.191) preserve the invariant measure with density /9 = — Dividing the second 

1-73 

and the third equations on the first and choosing a new independent variable, the angle of nutation 
6 = arccos73, we shall get the system of linear equations 



dKi _ J 3 sinfl K dK 2 _ mE 2 K x 

d0 s/hih + mR 2 ) 2 ' de ^h(h + mR 2 ) sin 6 

The general solution of these equations can be presented in the form [2 

I 3 sin 2 6> , + , n 1 — cos I 



(2.20) 



' 1 = c - ! vw^) f " +{ ' 1+ ^ 2 ' 

- ft - 7 =£*" ^1 + { ,i + „, *!±~f), P*) 

2V / i( / 3 + mR 2 ) L 

ts r< t?ic 1 1 — cos6K . ^ -. l + cos6> N 

K 2 = dF(£, r), 1, - ) + C 2 F{£, 77, 1, 



2 



where £ and 77 are the solutions of the quadratic equation x — x + = and F(£, 77, n, z) 

is the generalized hypergeometric function representable by series 



I 3 mR 2 
h{h+mR 2 ) 



x v 2 - r (£ + fe ) r 0? + k)T(n) z k 

Thus, the relations Q2.21JI define (implicitly) the integrals of motion. In this case they are the "con- 
stants" C\ and C 2 expressed through K\, K 2 , 9. 
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The quadrature for the angle of nutation can be obtained from the integral of energy written in 
the variables K\, K2, if 3, 9 



9 2 = 2 sin 2 6{h + mR 2 )P{9) 
P(9) = h 



K? 1 K? 



1 A 2 

— — mgRsmt 



(2.23) 



2/i sin 2 9 2 I x 

Here we assume that the variables K\ , if 2 are expressed through the constants of integrals and angle 9 
according to the formulas (|2.21|> . In this case the function P{9) (depending on the constants of 
integrals) define the analog of gyroscopic function for the Lagrange top [Hj, |2*2*j . 




Fig. 2. Phase portraits of the system l|2.2i}[l at various values C\ and C2. Left: the case of existence of three 
periodic solutions (C\ = 0.05, C2 = 0.01). Right: the case of existence of one periodic solution (Ci = 0.08, 
C 2 = -0.02). 

Thus, the equation l|2.23|) at the fixed values C\ and C2 define the one-degree Hamiltonian sys- 
tem. The phase portraits of this system on the plane 9, 9 are presented in fig. [2J All the variables 
73, K±, K2, if 3 are periodic functions of time with the period Tq and the corresponding frequency ujq. 

Remark 2. According to [1], the system (|2.19|) is the Hamiltonian one with degenerated Poisson bracket 
which has two Casimir functions expressed through hypcrgcomctric functions. 



2.4. Quadratures for angles of proper rotation and a precession 

According to the general Lie theory JS] , if the variables of reduced system (|2.18|l are the given functions 
of time, then all the variables of initial system (|2.12j) (and accordingly l|2.15jl ) can be obtained by one 
quadrature (if fields tfy (|2. 13|) and 1)2.14(1 are commuting). 

Indeed, using the equalities tg(/? = — (and correspondingly tg^ = — 7^) for angles ip and ip, we 
obtain 

73 Ki K 2 ; Ki 

* - — 3 t + Tfmrm - '* = (2 ' 24) 

Thus, for each of the angles the dependence on time is defined as an integral of a periodic function 
with the frequency ujg, hence it can be presented in the form (see, for example, JZJ; |24| ) 

(p = ui^t + (p*(t), tp = uj^t + tp*(t), (2.25) 

where <p*(t), ip*(t) are periodic function with frequency uog. Moreover, (|2.24|) and ()2.25|) imply also 
that all the frequencies ug, uj v , uj^ depend only on the constants of the first integrals. 
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Fig. 3. The surface of regular precessions. Parameters of the system are I\ = 0.25, h = 0.5, R = 1, m = 1, 
.9 = 1- 

2.5. A motion of the point of contact 

Following papers [2U], [21] we present the equation for the velocity of the point of contact in the form 

_^3_Ki _ K 2 _ \ e# (2>26) 

l- 7 3 2 / i y/hih + mR 2 )) 

where Z = X+iY and X, Y are the coordinates of the point of contact in the fixed frame of references. 

Thus the coordinates of the point of contact are determined by quadratures of quasiperiodic 
two- frequency (with the frequencies u>^, functions of time. 



Z = R 



3. The qualitative analysis and results 

Let's perform the qualitative analysis of the dynamics of the disk motion. We will make a classification 
of all possible motions depending on the constants of the first integrals. Some features of the considered 
case essentially complicate this work in comparison with the case of the Lagrange top for Euler-Poisson 
equations. For uniformity we recommend to study such analysis for the Lagrange case in book [3]. 
The complexity of analysis is caused by the facts that the integrals of motion can not be expressed in 
elementary functions (only in special one) and the system has no natural Hamiltonian presentation. 
Moreover, in addition to the motion of apexes of the body (disk) we shall classify trajectories of the 
point of contact obtained by additional quadratures of quasiperiodic functions. 

3.1. The bifurcation analysis of the reduced system 

Possible types of motion for the axis of symmetry of the body are completely determined by the form 
of the gyroscopic function P{6) and by the energy level. Critical values of the integrals of motion 
C\ , C2 , h are determined by the equations 

m = °' ~ir = °- (3 - 1} 

In three-dimensional space with coordinates C\ , C2 , h equations define a three-dimensional sur- 
face, so-called surface of regular precessions [S] (see fig. OJ). This name is connected to the fact that 
at the given values of integrals the coin performs motion with the fixed angle 9 = const, which is 
analogous to the precession for Lagrange top [22]. The full atlas of sections of the surface of regular 
precessions (bifurcation diagrams) by planes C\ + C2 = const and C\ — C2 = const is presented in 
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Fig. 4. Sections of the surface of regular precessions represented in fig. El by planes C\ + C2 = const. 
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Fig. 5. Sections of the surface of regular precessions represented in fig. by planes C\ — C2 = const. 

figs. 31 and [3 accordingly. In fig. and [7| for two different sections we show the forms of the gyroscopic 
function P(9) corresponding to various values of integrals C\, C2, h. 

Using these figures (and the rule of signs) we can easily study the stability of the corresponding 
solutions located on branches of the bifurcation diagram (branches corresponding to unstable solutions 
are represented on the diagram by a dotted line). In figs. El and vertical straight lines represent cases 
when C\ = or C2 = 0. In these cases the disk motion corresponds to the falling and planes 
determined by these equalities define in space of integrals C±, C2, h the two-dimensional manifold of 
fallings. Thus for almost all initial conditions the disk do not fall performing the rolling motion on a 
plane. 

Other remarkable motions correspond to the cases C\ = C2, the rolling motion of the disk, and 
C\ = —C2, the rotation of the disk about its axis passing through the diameter. During the latter 
motion the declination of the disk with respect to the vertical remains constant. 

Remark 3. The bifurcation diagram (fig. 01 El is different from one presented in papers [2], [231 since 
instead of the value of energy we use the value of angle of declination corresponding to the precession 0q and 
this function has no physical sense for other motions (when this angle is not preserved). Only the points on the 
surface of regular precessions have the physical sense. At the same time each value of constants C\, C%, h in 
space of integrals in fig. |31 corresponds to some motion whether this point is situated on the surface of regular 
precessions or not and this is important for the qualitative analysis. 

Remark 4. One of sections of the three-dimensional diagram by a plane h = const and the corresponding 
gyroscopic functions are presented in paper \7\. 
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Fig. 7. Various types of the gyroscopic function for the section of the surface of regular precessions by plane 
d - C 2 = 0.08. 

3.2. The qualitative analysis of motion of apexes 

The behavior of angles of proper rotation ip and precessions ip that together with 9 determined the 
motion of apexes is defined by relations (|2.25|) . The important feature in this problem is the two- 
frequency behavior of each of these angles. That is not usual for integrable systems. For example, 
for the Kovalevskaya top the angle ip{t) is defined by three frequencies jSj- In this case such phe- 
nomena is connected to the existence of two methods of reduction with respect to the symmetries of 
system (|TT2|) . (j2~T5|) . 

From the geometrical point of view whole space of variables M, a, (3, *y is foliated on three- 
dimensional tori defined as the joint level surfaces of the integrals C\, C2, h and the geometrical inte- 
grals. The motion represents a winding of the three-dimensional torus with frequencies ujq, uj v , uj^ |24j . 
(For the reduced systems (j2.12jl and (|2.15jl the corresponding tori are two-dimensional.) 

Since the frequencies depend only on constants of the first integrals, all motion on the torus 
have the identical frequency that not evident for nonholonomic systems. Even for the integrable 
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Fig. 8. Trajectories of the point of contact of the disk in the absolute space at various values of the integral 
of energy. Parameters of the system correspond to figure a). The closed trajectories in the frame of references 
rotating with the angular velocity uj^ (see explanations in the text) are presented in the right upper corner of 
each figure (except the case of infinite motions). In figures a) and b) we present various types of motion of the 
disk at the energy h = 0.86. Figures c) and d) correspond to the energy h = 0.92217 when one of the motions 
becomes resonant (uig = uj^ ) and the secular drift (fig. d) is observed. The increase of the energy in figures e) 
and f) to h = 0.961 makes both types of motion bounded again. In figure g) the motion of the disk is presented 
at h = 1.1 after merging of two domains of possible motions corresponding to various types of motion. The 
infinite motion in figure h) corresponds to the resonance uj^ = 2ue at the energy h = 1.18169. In figure i) the 
motion of the point of contact of the disk is presented after the further increase of the energy up to h = 1.4. 



nonholonomic systems on two-dimensional tori there is a non-uniform rectilinear motion and, generally 
speaking, the intermixing is possible (see paper [1]). 

Practically these arguments prove that the given system is Hamiltonian one in the analytical sense 
( though the Hamiltonian function can be a different from the energy (|2.9ft |24j . Moreover, taking into 
account only the analytical point of view we can say that near the nonsingular torus the system becomes 
the Hamiltonian one by the infinite number of methods JS] • 

Remark 5. N. K. Moschuk in (2H1 observed a related phenomena studying the nonholonomic Chaplygin 
system possessing some number of the linear with respect to velocities first integrals. 
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Fig. 9. Dependencies of frequencies mg, and w v from the energy at I\ = -r, I3 = ^, R = 1, m = 1, 
g = l. and various values of integrals Ci, C2. In figures a) and c) the areas marked by the rectangles are 
separately presented in the increased scale. In the field of energies where two different types of motions 
are possible in the absolute space we denote the frequencies by w^,, and w^,, u^. The resonance 
energies are marked on the graphs by thick dots. The orders of the resonances are indicated near the 
dots. The values of integrals for the dependencies presented here are: a) C\ = 0.04, C2 = —0.02; 
b) Ci = 0.09, C 2 = -0.07; c)d = 0.065, C 2 = 0.055; d) C x = 0.09, C 2 = 0.03. 



At the same time the existence of a natural (algebraic) Poisson structure with a Hamiltonian 
defined by the energy (|2.9|) remains an open problem. A. V. Borisov and I. S. Mamaev show that the 
reduced system ()2.19j) is the Hamiltonian one with some algebraic nonlinear bracket (see (5]), however 
the possibility of its lifting on the systems (|2.12|) and 1)2.15(1 is still not investigated. 



3.3. The analysis of motion of the point of contact. 

For the analysis of motion of the point of contact we decompose the velocity (|2.26|) in the Fourier 
series with respect to time. Then from (j2.25|) we get 



Integrating with respect to time we obtain 

Z(t) = Z + e^* 



-* i{u^ + TWJq) 



Thus, if at + nujg / Owe use the frame of references rotating about the point Zq with the angular 
velocity then the point of contact traces some closed curve (see (21], (20] )• Various types of such 
closed curves and trajectories corresponding to them in the fixed space are presented in figure [S] 

At the resonance aty + uluq = we observe the secular drift of the point of contact. Graphs of 
frequencies co^(h), co$(h), u v (h) at the fixed values of integrals Ci, C2 are presented in fig. O They 
show that the relation uj^, + ncog = can be fulfilled both in the case of existence of one and of three 
regular precessions. And at the same energy some initial conditions lead to a secular drift while the 
others are not (see fig. |H]). Since all frequencies depend only on the values of the first integrals the 
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relation uty + uujq = define in three-dimensional space of integrals some two-dimensional manifold 
corresponding to the infinite trajectories of the disk. 

Thus, for almost all initial conditions {except the indicated manifold) all trajectories of the disk 
are bounded. 

We can consider this result to be opposite to the one obtained from research of dynamics of 
the point of contact for the Chaplygin ball on a horizontal plane (see ^2]) where the majority of 
trajectories, on the contrary, were unbounded. 
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